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O ■ We argue that recent experiments by Kirtley et al. pl may show evidence of 



time reversal symmetry breaking in YBa2Cu307 at crystal grain boundaries. 



C^ ■ We illustrate this through a Ginzburg-Landau model calculation. Further 

experimental tests are proposed. 
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In a recent paper, Kirtley et al. |jT[| reported the observation of magnetic defects at ar- 
tificially engineered grain boundaries in thin films of the high temperature superconductor 
YBa2Cu307 (YBCO). The grain boundaries were the borders between a triangular YBCO 
inclusion in a film of YBCO with the crystal axes misoriented with respect to one another 
in the two domains (inside and outside of the triangle). While the resolution of the mag- 
netic microscope (~ lO/im which is roughly ten times the estimated Josephson penetration 
depth A J 0) used for detection is not sufficient to tell with absolute certainty, the observed 
magnetic defects appear from their shape and localization to be superconducting vortices 
carrying small fractions of a flux quantum $o = hc/2e. These vortices are attached mainly 
to the corners of the triangle, but occasionally appear along the edges of the triangular 
inclusion. The purpose of this Letter is to point out that the identification of these defects 
as fractional vortices, if correct, demonstrates that the materials in question have super- 
conducting order parameters, and thus ground states, that violate time reversal symmetry. 
The experiments cannot tell whether the T-violation is a bulk or interface (grain bound- 
ary) effect. Our argument is simply that the flux carried by a vortex measures the phase 
defect of the order parameter along a closed path encircling the vortex, and can therefore 
be fractional only if the order parameter undergoes a phase jump A0 not a multiple of 2% 
along this path. The recent Josephson tunneling experiment of Wollman et al. [^ and the 
observation of half-integer flux quanta by Tsuei et al. [Q are speciflc examples of this for 
which A0 is an odd multiple of vr. Because of speciflc symmetry properties of the Josephson 
junctions their results were interpreted as strong evidence for (ij.2_j^2-wave pairing symmetry 
in YBCO, which is a T-conserving superconducting state. On the other hand, the recent 
experiment of Kirtley and co-workers [|l| can be explained only if A0 is not a multiple of vr, 
which requires T-violation. 

Let us briefly review the historical context of T-violation in unconventional (both heavy- 
fermion and high-Tc) superconductivity. It has long been suspected that time reversal sym- 
metry breaking is responsible for some of the unusual magnetic properties of heavy fermion 
superconductors, in particular (U,Th)Bei3 and UPts [|],0. The possible appearance of frac- 



tional vortices in these materials has aheady been suggested [§ and investigated theoreti- 
cally ip. The conditions under which a superconductor with a real order parameter in the 
bulk phase may spontaneously break time reversal symmetry have also been studied in the 
context of Ginzburg-Landau theory |TO| , p!T| . Surfaces and domain walls were found under 



certain conditions to favor the formation of a locally T-violating state as a means of low- 
ering the energy cost of an inhomogeneous order parameter [|ll|]. T- violation (specifically, 



a (ia;2-j/2 + iedxy order parameter) has been predicted in high- Tc superconductivity via the 
anyon technique applied to the t-J model [0,|ll]- A superconducting state with s + idxy 
symmetry has also been proposed [l^. However, none of the telltale signs of T- violation has 



been detected in bulk measurements [jT5[. This does not preclude the existence of a complex 
order parameter at surfaces and grain boundaries since bulk measurements are not sensitive 
to the existence of such a phase. It is an additional point of this Letter to show at least on 
a phenomenological level that such states are indeed possible. 

To illustrate our idea we first analyze the properties of superconducting states near an 
interface by means of a Ginzburg-Landau (GL) theory of two complex order parameters, 
?7i and 772. These two order parameters belong to pairing states of different symmetry, e.g., 
dx2-y2 and dxy'- V'i(k) = k^ — k^ and 4'2(X) = k^ky, which are non-degenerate under the 
tetragonal (-D4/1) crystal field symmetry assumed here. From this we can derive a GL free 
energy functional of rji and 7^2 with the requirement that it be a scalar under all symmetries 
of the system (for a review see Ref. 0): F = Fi + F2 + F12 with 

F, = J d'x [«.(T)|r/,|2 + A|r/,|4 + K,\Bv^\^] (1) 

F12 = / d'x b\vi\'\V2\' + Sivfvl + Vfvf)] (2) 

where aj(T) oc T — T^ {T^, the bare bulk transition temperature of the order parameter 
rji) and (3i,Ki,'^ and 5 are real phenomenological parameters which contain all the relevant 
physical information of microscopic origin. The gradient terms are given in the gauge in- 
variant form D = V — i27rA/$o with A the vector potential (B = V x A) and $0 the flux 



quantum hc/2e. We note that this is not the most general Ginzburg-Landau free energy 
allowed by symmetry. For simplicity we include only terms which are relevant for our dis- 
cussion and, in particular, use an isotropic gradient term. The boundary conditions at the 
interface (surface or grain boundary) can be formulated in the standard way by 

n ■ D?7, = r]i/hi, i = 1,2 (3) 

at the interface. Here n denotes the normal vector of the interface and bi is the so-called 
extrapolation length depending on the properties and orientation of the interface (see PJ^). 
Non-s-wave order parameters are often distorted in this manner in the vicinity of an interface 
due to scattering effects p. 

Let us now examine the possibility of T-violation at the interface. We shall assume, as 
suggested from various experimental observations, that in the bulk only the single component 
rji = uic^'^^ exists, while 772 = U2e^'^^ vanishes for all temperatures. Thus, we require that 



Ul 



H = u{T) = 7-ai/2A, U2 = (4) 

which is satisfied for all temperatures below T^ under the conditions 

(7 + 2(5cos(2^))m2 - a2 > 

(5) 
Tc, > Tc, 

where u is the asymptotic value of Ui in the bulk region and = 4'i — 4'2 denotes the relative 
phase between the two order parameter components. With the choice 6 > 0, the state with 
6 = ±7r/2 (ipCk) = uipi{\i.) ± if-?/'2(k, or dx2__y2 ± ied^y) is closest to the instability although 
under condition, Eq.(5), not stable for any temperature. Thus, we may first treat 'Ui(x) near 
a planar interface as though U2 were zero, considering the GL differential equation obtained 
by variation of F with respect to rji (we neglect the vector potential for this discussion). 
Note that ui depends only on the coordinate parallel to the normal vector n which we may 
choose to be parallel to the x-axis with the interface located at x = 0. 



Ki-g^ = a,u, + 2(3,ul (6) 

As a solution to this equation we obtain 

Ui{x) = -utanh — - — (7) 

with ^ = J2Ki/ai and xq = (^/2) sinh~ (46i/^). Next, we ask whether this interface state 
could be unstable against the admixture of a small component U2- This question can be 
answered by analyzing the linearized GL-equation of U2 for fixed ui{x). 



d'^U2 

dx"^ 



K2^ = a2U2 + (7 + 2(5 cos{29))ul{x)u2 (8) 



It is easy to see that this equation has the form of a Schrodinger equation for the wave 

function U2 of a particle in a potential well (including the boundary condition for U2). The 

"lowest energy eigenstate", which is a bound state, defines the critical temperature T* 

below which the interface state, Eq.(6), is unstable. The corresponding wave function U2{x) 

is nodeless and decays exponentially in the bulk region. For 5 > 0, as assumed above, the 

relative phase 6 is ±7r/2 so that this state breaks time reversal symmetry and is two-fold 

degenerate. It is not possible to obtain an analytic solution of Eq.(8) in general. However, 

under the rather restrictive condition 62 = ~^coth(xo/^) the instability condition can be 

given analytically as a2{T*) + m^(T*)(7 — 26) /2 = with the bound state wave function 

const. . , 

"^(") = cosh((x + xo)/0- ^'^ 

Our analysis demonstrates that under certain conditions the interface of an unconventional 
superconductor can give rise to a locally T- violating state (see also |jll|). Furthermore, T* 
(< Tci) defined above is the temperature at which a continuous phase transition from a 
T-conserving (T > T*) to a T-violating state (T < T*) occurs. It is, however, not our aim 
to discuss here a possible microscopic basis for our GL-model. Rather we are interested in 
some of the consequences of a T-violating superconducting phase. 

Let us now study the phenomena which occur at a Josephson junction between two su- 
perconductors A and B if T- violation is present. The following discussion does not depend 



on whether the T-violation is a bulk or, as discussed above, an interface (junction) phe- 
nomenon. Because we have two complex order parameters at the interface, the Josephson 
phase-current relation consists of four terms 

J = J2 Jij sin(0jB - 0ja) (10) 

where Jij are real constants whose sign and magnitude depend on the grain orientation 
and order parameter magnitude at the interface: Jij oc \rjiB\\rjjA\Xiip-B)Xjip-B), '^a,b is the 
junction normal vector on either side and, typically, Xi(n) = "^x """-^ ^ind X2(ii) = nxHy. We 
assume that the current through the interface vanishes, because due to screening effects (on 
a length scale Aj), such currents can only flow near the boundary of the interface or near 
a vortex. Furthermore, we assume that the couplings Jij are sufficiently weak so that the 
relative phase between rji and rj2 is not affected, i.e.: (piA — (p2A = 'Pie — 02B = ±vr/2 in the 
T-violating state. This simplification is not important for any of our later conclusions and 
a more complete discussion will be given elsewhere. 

The latter assumption allows us to minimize the junction energy, E = 
— ($o/27rc) J2i,j Jij cos(0jB ~ (pjA), by choosing the phases such that J = 0. We obtain 

A0, = 0,B - 0.A = ± tan-1 f ^^^— ^) (11) 



for a junction a with all Jj,- > 0, and 



Ji2 — J21 



^4>b = (t>iB - 0jA = TT ± tan M - — --— (12) 

V Jii + J22/ 

for a junction h with J12, J12 > and J21, J22 < 0. 

We consider now the situation where these two types of junctions, a and 6, intersect 
(forming a grain boundary corner). Such a corner is accompanied with phase winding or 
a vortex, because, in general, A^^ 7^ A05. For the calculation of the magnetic fiux of this 
vortex we notice that the supercurrent is given by the expression 

^j= E^^«^(Vv^-|^A + V0.) (13) 

C i=l,2 *0 



with Tjj = Uje^^'^^^'^\ j=l,2, and y? a phase of the order parameter continuous even at the 
junction and < 0j < 27r . We choose a path C encirchng the corner at a distance far 
enough so that j = on C. We denote the segments of C in superconductors A and B by 
Ca and Cb respectively. Using 0i = 02 ± 7r/2 the circular integral of Eq.(12) on C leads to 
the flux 



^=n+{f ds+ f ds) ■ V01 = n + ^^\ "^" (14) 

$0 Jca Jcb 2% 

where n is the integer winding number of if. Obviously, the flux at the corner can have any 
fraction of $o and is determined only by the properties of the junctions. On the other hand, 
it is easy to see from our discussion that in the case of a T-conserving superconducting state 
the only fractional vortex is that with half a flux quantum $o (*^ = *^o('^ + 1/2)) Ill6|,p|. 
The field distribution of such vortices would extend along the junction on a length scale Xj 
while penetrating the bulk only by the London penetration depth A <^ \j. 

The twofold degeneracy of the T- violating interface state implies the existence of domains 
and domain walls. There is a phase winding and flux associated with the intersection of a 
domain wall and a grain boundary, because the phase jump A0 at the junction is different 
on the right and left hand side of a domain wall. Following above scheme, a domain wall on 
junction a contains a flux $/$o = n + A^q/tt. These vortices are similar to the fractional 
domain wall vortices analyzed in Ref. 0. They are not connected with corners, but can 
essentially by located anywhere on a grain boundary. Hence, we may conclude that our 
model can account for both fractional vortices at the corners and along the edges of the 
triangle as observed by Kirtley et al. |]l[. 

Let us make several remarks. The central point of this work is that any superconducting 
state with fractional vortices containing other than (n + 1/2) flux quanta violates time 
reversal symmetry. Fractional vortices are not specific to the d^2_y2 + ied^y state, although 
in the interest of simplicity we restricted ourselves to this order parameter in our model 
calculations. It should also be noted that on the basis of the experiment by Kirtley and 
co-workers alone, the specific form of T-violation cannot be deduced. Our conclusion that 



time reversal symmetry breaking has been seen can only be wrong if the assumption of 
existence of the fractional vortices is wrong. There are several experiments which would add 
considerably to our understanding of this matter. One such experiment would be to look for 
the critical temperature T* at which T-violation occurs. According to our discussion there 
would have to be a second phase transition below the onset of superconductivity, although 
this could very likely be a grain boundary phenomenon only. Above this temperature T* 
there can be no fractional vortices apart from those with $ = ±$o/2. It would also be 
interesting to look for fractional vortices in different materials such as HgBa2Ca2Cu308 or 
Bi2Sr2CaCu208, perhaps utilizing different geometries which might better isolate the grain 
boundary corners. The existence of a complex order parameter at the interface may also 
shed new light on the interpretation of Josephson junction experiments such as those of 
Chaudhari et al. [|1^] and Sun et al. ||18[, beyond the analysis given recently by Millis 0. 



In this Letter we have shown that (1) the existence of vortices enclosing a fraction of a 
flux quantum requires the breaking of time reversal symmetry, and (2) that the converse is 
also true. We argue that this has been observed at grain boundaries in YBa2Cu307. Further 
experiments are needed to deduce the nature and extent of T-violation in high temperature 
superconductors. 
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